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The one dimensional Dirac equation with a rational potential is reducible to a ordinary differential 
equation with a Riccati-like coefficient. Its integrability can be studied with the help of differential 
04 . Galois theory, although the results have to be stated with recursive relations, because in general 

the equation is of Heun type. The inverse problem of finding integrable rational potentials based on 

■ the properties of the equation's singular points is also considered and two classes of potentials are 

<n: recovered. 

< 

<n: 

■ The goal of this work is to examine the solvability of the Dirac equation, when reduced to one spatial dimension. 

(-H ' The key point is of course that the analysis is easier and deeper for an ordinary differential equation, and also that the 
O J reduction itself can be achieved by many means, such as separating the radial from the angular part, or neglecting two 
directions in the Cartesian coordinates. A natural source of the potential is the magnetic field, which also provides a 
direct connection with systems of quantum optics and a possible experimental realization [l| . The resulting equation 
Ci • form is also encountered when dealing with the supersymmetric Shrodingcr equation [3], to which all results on 
^ ] integrability can be transferred. Of course, one has to keep track of the different definitions of potentials and physical 
interpretation, but the mathematical structure remains the same, and that is the main interest of the present paper. 

A great many of known solutions were obtained for the scale invariant potentials, and can be found in Q. However, 
such approach relies on the assumption of a particular symmetry and is not exhaustive. Not disregarding it, it seems 

■ natural to also consider equations like this (i.e. second order and linear) on grounds of the differential Galois group 
T^lj- , theory. Mostly because it provides rigorous tools for finding the exact forms of solutions and also because it strictly 
T— I ■ rules out unsolvable cases - exactly like the "ordinary" Galois theory for polynomials. 

'', Since linear equations often present themselves as central, like the Schrodinger equation or the density perturbation 

■ in cosmology, they have already been successfully studied with these methods, e.g. in 0, S Hi- When it comes to 
[ the Dirac equation, the complete analysis in the case of general polynomial potentials was carried out in a previous 
. work by the author Q . This class happened to be very restrictive, and only one non-trivial case - the so called Dirac 

^""^ ■ oscillator - is explicitly integrable. 

^ . The next obvious step seems to be the wider class of rational functions which is undertaken here. This turns out to 
be too wide to be treated as completely as polynomials, the basic reason being that it is the poles of the equation's 
coefficients that matter, and polynomials could be said to have just one - at infinity. As will be shown, even the simplest 
^ ■ potential with just one pole and linear polynomial part cannot be analyzed explicitly in that it requires recursion to 
■ ■ ' formulate integrability conditions. This is a typical feature of the Heun equation and its generalizations, and thus a 
generic rational potential is beyond explicit classification which was possible for polynomial ones. Fortunately, the 
recursive conditions can be studied to some extent with regard to the real values of the parameters, and shed some 
light on the number of admissible physical cases. 

The analogy with polynomials also motivates the study of potentials which give rise to similar behaviour of the 
equation's singular points and since direct analysis is impossible, one might also try the inverse problem. Thus, 
additionally, two other scenarios are considered here: rational potentials which produce only the irregular singularity 
at infinity, and the Whittaker equation, whose solvability is known and can be used to obtain a new, two-parameter 
family of integrable potentials. 
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II. 

To be more specific, let us consider the same setup as in ^], so that the time-independent Dirac equation in one 
dimension z is 

E^^{-iazd,^^Tn)i,, (1) 
which gives the second order equations for each spinorial component 

^'IM = [y{^f ± y\^) - E^) ^lA^), (2) 

where the Dirac matrices were taken to be 

and the potential V was introduced in the so called scalar coupling, so that it is added to the mass term m ^ m + V . 
For simplicity m will henceforth denote the effective mass thus obtained, that is the old mass plus the constant term 
of V . The other, vector, coupling is equivalent, as shown in iS], so that essentially the equation in question is just 

'4," = {V^ + V' - E'^)ij =: r i). (4) 

The question of exact solutions of such an equation with a rational coefficient r{z) G C(z) is best expressed in terms 
of field extensions of C(z). An extension including all the elementary functions, direct integration and some special 
transcendental functions is the generalized Liouvillian extension. It consists of a finite chain of differential function 
fields, 

FoCFiC...CF„, (5) 

such that each new field Fk+i is obtained from the previous Fk by adjoining an element Xk which satisfies one of the 
three following conditions: 

1. Xfe is algebraic over Fj,. 

2. Xk is primitive, i.e. x'j, £ Fk- 

3. Xk is the solution of x' ~ ax, for some a in Fk- 

These fields are differential fields precisely so that there is some notion of the derivative x'j. , although here it is simply 
the differentiation with regard to z and the starting field Fq is just C(2;). If the solution of a given equation belongs to 
Fn we might call it a Liouvillian solution and it is explicitly expressible by finite number of elementary operations like 
integrals, exponentials or solving algebraic equations. In other words, the equation in question is solvable/integrable 
if it has Liouvillian solutions. 

A more detailed exposition of the above can be found in , and here it should be added that this particular notion 
of exact solution is intimately connected with the underlying group theory. Namely, one defines the differential Galois 
group of a differential equation as the group of automorphisms of the field to which the solutions belong, such that 
their action on the base field Fq is trivial. Intuitively it means that the group permutes the solutions, but does not 
change the equation's coefficients. It then follows that the equation is solvable when the group is solvable. In case 
of the second order equations like ^ the group is an algebraic subgroup of SL(2,C), and there exists an algorithm 
contrived by Kovacic [8] for deciding on its integrability and finding the solutions. Also, most of the common equations 
of mathematical physics have been completely analyzed in this regard (see e.g. the monograph for a collection of 
results), so the question is often that of finding a connection with what is already known. This will be the main route 
taken here, allowing to uncover a new, wide class of integrable potentials. 

III. 

A naive generalization of the results of @ is difficult for the following reason. When V is polynomial, it only has 
one pole at infinity (and so has r{z)), and it just the Laurent expansion with finitely many terms. The polynomial's 
coefficients are then all we need for the analysis. A rational V can be defined by a quotient of two polynomials, and 
although the denominator still has finitely many coefficients, they are not used directly. It is the poles that matter 
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so r{z) needs to be expanded around each one. One is therefore left with an arbitrary number of poles of arbitrary 
orders and coefficients with which to construct a general rational potential that could be analysed effectively with 
Kovacic's algorithm. 

Unfortunately, even the seemingly simplest case of adding just one pole to an integrable polynomial potential leads 
to only implicit integrability relations, so a full classification analogous to the polynomial case is impossible. To see 
this let us consider the following 

V{z) = - +m + Xz, (6) 

with a ^ 0. In general this yields the (biconfluent) Heun differential equation whose solvability is determined by 
algebraic relations defined only recursively - the subject of this section. Explicit analysis is possible whenever one of 
the other parameters is zero, and the equation reduces to a special case of the confluent hypergeometric one, which 
has countable families of Liouvillian solutions {oj . 

More specifically, for A = ([5]) becomes a special case of the Whittaker equation, analysed in detail in the next 
section. When A 7^ but to = it reduces to the Laguerre equation 

1 

ug"{u) + (a + - - u)g'{u) = (7) 

where 

z-°V;(z) = e-"/2g(-y)^ w = -Az^. (8) 

The solutions are the Laguerre functions ^^'^\u) with n = and these are again expressible by confluent 

hypergeometric functions. They are Liouvillian if at least one of n or n + a — i is an integer (with one solution 
polynomial for positive integer n). Additionally, E = always leads to an explicit solution of 

V' = exp J Vdz, (9) 

regardless of the form of the potential. 

The generic case with all parameters non-zero can be slightly simplified with taking A = ±1 or, in other words, 
rescaling as follows 



z^z/V|A|, to->toV|A|, £:^£;v|A|. (10) 

The reason for using absolute value is to keep the physical quantities real. This will make identifying the physical 
integrable cases easier. 

Turning now to Kovacic's algorithm, the first case will be the only one to possible hold, and the solution must be 
of the form 



ip = P{z)exp J uj{z)dz, (11) 

with rational w and polynomial P. Denoting ip' /ip —: and substituting into ([5]) we have 

e' + 0"^ = V^ + V' - E'^, (12) 

so that the poles of V and 9 must agree except that there could be additional simple poles of 9 with residue 1 not 
appearing in V. Since 9 = -p- -t-w, the first term is needed precisely to account for those poles, and w has to have the 
same polynomial part as V (up to a sign) so that the expansions at infinity agree. Finally, the pole a/ z requires that 
uj includes either a/ z or (1 — a)/z, so that there are four possibilities all in all: 

a 1 — a 
wi^2 = — ± (m + Az), W3_4 — q= (to + Xz). (13) 



The next step is to substitute the above into the equation in question and look for polynomial solutions for P with 
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its degree d determined by the highest order term in each respective equation. The four subcases are thus 
1. P" + 2(-+m + Xz)p' = 2XdP, d = - — 



P 



(14) 



3. p" + 2 ( ^ -m-Xz\p' ^(—- 2\d\ P, d=^-2, 



z 



4. P" + 2i'—^+m + Xz] P' ^ 2X^ R d=-£ + 2«-l. 



2A 

E 
2X 



Taking 



P = Y.Piz\ (15) 



1=0 

the above can be stated as recurrence relations 

1. l{l + 2a - l)pi + 2m{l - l)p;-i = 2A(d - / + 2)p/_2, 

2. l{l + 2a - l)pi - 2m{l + 2a - l)p;-i = 2A(/ -d- 2)p/_2, 

3. + 2a)pi - 2mlpi^i = 2X{1 - d - 2)p/_2, 

4. Z(; + 1 - 2a)j5/ + 2m{l - 2a)p,_i = 2A(d - I + 2)p,_2, 



(16) 



with common initial conditions of 



Pd+i =0, pd = l- (17) 



They reflect the fact that P can always be taken to be a monic polynomial, and its degree is exactly d. Consequently, 
the relations hold for ^ < d + 1 and one can successively generate p;, as functions of the parameters, down to pq. The 
last but one relation, i.e. for I — 2 gives us po but there is one more with I = 1 and only two terms pi and po since for 
polynomials p-i = 0. This is the required condition for the existence of polynomial P. It is convenient, however, to 
consider the last relation as the definition of which will be a function of to. A, a and d, and then the condition 
p_i = (equivalent to the two-term relation between pi and po) will yield the admissible values of parameters for 
which a polynomial P exists, and hence the main equation is solvable. Since the whole problem is considered over the 
complex domain, all polynomials have the maximal number of roots, counting with multiplicities. When considered 
as polynomials in m, p_i is of degree d + 1 provided that the middle term of the recurrence does not vanish. This 
particular way that to enters also means that all the polynomials are either even or odd functions of m. 

For each degree d the relations (|16p implicitly define p_i whose roots then provide the value of parameters leading to 
an integrable case, while the energy is determined by its relation with d in each case. Of course, one is mostly interested 
in the real values of the parameters as they correspond to the physical case. This is the reason for approaching pi as 
polynomials in to and a (A can be scaled away, as mentioned before). The particular form of the recurrence makes 
those sequences define orthogonal polynomials in most cases (although this property is of no use here), and their zeros 
are then all real, simple and with the separation property 

(k) 

where ai is the z-th zero of the A:-th polynomial. To see this in the first subcase note that 

Pd-i = -dm, af = 0, (19) 

and 



p._2 = d[2[d - IW + 1 - 2a - d)/4, ar^' = -f-^^^ -'t'' = \/ (^0) 

so that distinct real roots exist for 2a > 1 — d. c? = is ruled out by ii^ 7^ 0, d = 1 ruled out by m ^ 0, so the smallest 
d is 2, which translates to a > —1/2. Assuming that the separation property holds for polynomials pi and p/_i the 
proof is by induction on I. Evaluating pi-2 at a zero of gives 

l{l + 2a - l)pi{a^t^^) = -2(d - I + 2)pz_2(af "'^), (21) 
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where A was taken negative, as is necessary by the energy condition in (|14l) . Since the zeros of pi lie only between 
those of pi^i it fohows that the left hand side of the above is non-zero and switches sign as i increases, provided that 
Z + 2a — 1 is non-zero. It follows that pi-2 must change sign in each interval (af ^\ ^i+i^''*) which can happen only at 

a root. The only thing left to show that p/_2 has two roots outside of (of a|j_j^-^). With the further assumption 
that Z -I- 2q; — 1 > the sign of p/_2 at the last root of pi-i is opposite to the sign of pi but they both have the same 
sign in their leading powers of m which are m'*"''^ and m'*~^ respectively, so that pi-2 must have another root for 
the signs to agree as m — oo. The same happens at the lowest root. Finally, care has to be taken as to the middle 
term of the recurrence, as it cannot vanish (the number of roots must be maximal). Fortunately for I — 1 the relation 
just states that api = 2X(d + so effectively we can take / > 3 and pi = as the integrability condition, and 

automatically the smallest d to consider is 3 as d = 2 implies p_i = —2m. Altogether this means that the induction 
is valid for a > —1, and the separation property holds up to pi. 

The third subcase is analogous and even simpler, with the only difference being that the energy condition requires 
A > 0, and the positivity of ^ -I- 2a — 1 is enough to guarantee the appropriate changes of signs, and maximal degree 
of each pi. For the first induction step, the condition for pd-2 to have simple real roots is that d — 2a -I- 1 be positive. 
Since d — gives p^i — m, it follows that a > 1 is enough for the separation property. 

Subcases two and four differ in that the sign of A is now arbitrary, but the general proofs are the same and we get 
A>0, a>OorA<0, a< —d/2 in subcase two, and A > 0, a > (d -I- l)/2 or A < 0, a < 1/2. Where again the 
trivial cases to = were discarded. However, the results can be slightly strengthened here. In subcase two, for a — 
the relation with I = 1 just gives p-i = 0, and one can see that for — 1/2 < a < 0, the induction can still be carried 
as before, with the modification that signs at the last step are reverted. Namely 

2api - 4TOapo = -2A(o?+ l)p_i, (22) 

so that for A > at the maximal root of pq the signs of p_i and pi are the same, but additionally the leading coefficient 
differ in sign, as the coefficient multiplying mpQ in the relation is positive, while it was negative for all higher I. Thus 
the existence of the highest root follows. For A < this cannot be done, as the main restriction comes from the roots 
of pd-2 not from the third induction step. 

The same is true for subcase four with (d + l)/2 > a > d/2 and A > so that the region where the separation 
property holds can be extended to a > d/2. The case of a = (d + l)/2 is, by the energy condition, equivalent to 
-E = so it brings no new solutions. 

Putting the above together, and discarding the m = case we have 

Theorem 1 The integrability condition p-i =0, considered as an equation in m has the full number of single real 
roots in the subcases (|16p for 

1. X <0 and a > —1. There are [(d— 1)/2J positive roots with d ~ 

2. X > and a>— -^j or X < and a < —d/2. There are [d/2] positive roots with d — — 2a — 1. 
5. A > and a < 1. There are [d/2] positive roots with d — 

4. X > and a > d/2, or X < and a < ^. There are [d/2] positive roots with d = + 2a — 1. 

Next, the inspection of the to = case tells us that for odd d = 2n + 1 the conditions read 

^^fr'^ (2n+l)(a + n) 
P-'= 11 A(d-2n-hl) 

J^t^^ (2n+l)(a + n) _ 
li A(2n-d + 3) 

n—0 

-li A(2n-d + 3) 

n—0 

^^f.^^ (2n + l)(n + l-a) 
11 A(d-2n+l) 

n—0 



(23) 
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so that there are always n = (d — l)/2 real simple roots. This happens because the recurrence involves only two 
terms then, and we start with reducing the index by 2 at each iteration to get to It then follows that in some 
neighbourhood of m = there must still be n real simple roots, because p-i is a polynomial in a whose coefficients 
are again polynomials of to, so that the roots must be continuous functions of m. Note however, that in subcases 2 
and 4, the full polynomial could be of higher degree in a, thanks to the middle term, and there will be additional 
roots, that go to infinity when to — > 0. 

Even degrees are much more involved, so we include the proof in the appendix, and state the general result here 

Theorem 2 In the vicinity of m — there is always the maximal number of real simple roots in a of the integrahility 
condition p_i = 0. Depending on the parity of d and the subcases of (jl6p one has 

1. For odd d there are [d + l)/2 roots (including a — 0) given by ([23]). For even d there are d/2 roots interleaved 
with the roots of (jA3p . 

2. For odd d there are (d + l)/2 roots (including a = Q) given by (|23p . For even d there are d/2 roots interleaved 
with the roots of (|A19P ( also including a = 0). 

3. For odd d there are (d + l)/2 roots given by (j23l) . For even d there are d/2 roots interleaved with the roots of 

(Em. 

4-. For odd d there are (d + l)/2 roots, one of them being a = (d + l)/2, given by (|23p . For even d there are d/2 
roots interleaved with the roots of (jA22[) . except for the the shared root a = {d + l)/2. 



The the best way to view all the branches of the algebraic sets defined by the integrahility conditions, and to get 
the intuition about the asymptotic behaviour of roots, is a plot of p_i = in the (to, a) plane. 




-4 -2 2 4 -3 -2 -1 1 2 3 



m m 

Figure 1: The curves (solid) defined by pi = in subcases 1 and p_i = in subcase 3, both with d — 7. Circles represent 
zeroes obtained with the separation property, while squares are zeroes for m = that change into simple roots for small positive 
(and negative) m. 



The graphs suggest that there are horizontal asymptotes in subcases 2 and 4, so that for large to there would also 
need to exist a fixed number of solutions. To check this, let us take rescaled polynomials defined by pi — m'^'^^ri which 
the satisfy 



2. l{l + 2a- l)-^ - 2{l + 2a- l)r;_i = 2A(/ -d- 2)r/_2, 
4. l{l + l- 2a)-^ + 2{l - 2Q)n_i = 2X{d - I + 2)r,_2, 



(24) 



where the rescaling is based on the maximal power of to appearing in pi . This means that each r; is finite as m — > cx) 
and the limit case of the relations reads 

2. -2(/ + 2a-l)n_i =2A(?-d-2)n_2, 
4. 2(?-2a)n_i =2A(d-? + 2)n_2, 
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so that at each recurrence step a real zero in a is introduced - exactly as in the m = case above. In subcase 2 
I = 1 adds the precluded a = 0, while in subcase 4, Z = d + 1 is the previously mentioned straight line of solutions 
(independent of m). These real roots at infinity cannot turn into complex one for large finite m since any complex 
roots would have to come in conjugate pairs. This would require the polynomials of twice the degree, but also the 
roots at infinity to be double. By continuity we thus obtain 

Theorem 3 For sufficiently large there are always solutions o/p_i = 0, when considered as an equation in a. In 
subcase 2, there are d solutions which asymptotically approach a = {1 — l)/2 for I = 2, . . . , d + I. In subcase 4, there 
are d solutions that approach a = 1/2 for I = 1, . . . ,d and the line a — {d + l)/2 is also always a solution. 

The question remains as to the vertical asymptotics of the branches, that is what happens for infinite a. Rescaling 
this time by pi = a'^~Wi gives 

2. -4TOn_i = 2A(Z-d-2)n_2, 
4. -4TOn_i = 2A(d-Z + 2)n_2, 

so that only the factor of m is introduced at each step. This means that all the branches tend to m = 0. For any 
straight line of constant m, each branch must intersect it then, as the previous theorem guarantees that they extend 
towards infinite to. If in the vicinity of to = it were possible to prove the existence of solutions also for odd degrees, 
that would translate into global existence of exactly d + 1 roots in a for any value of to (counting also a = and 
a = (d + l)/2 where necessary). 

Finally, the last two cases of Kovacic's algorithm need to be dealt with. The second case cannot hold for the required 
quantities Ec as defined in [8| are: Eq — {4(1 — a), 2, 4a} and Eoo = {—2}. Then the degree of the polynomial P 
defined as above should be d = (coo — &o)/2 for Cx G E^, and this implies that a has to be an integer or half integer. 
However this makes all Cx even, so that case two cannot hold. Case three is ruled out because the order of r{z) at 
infinity is less than 2. In other words infinity is an irregular singular point. 



IV. 



The direct approach above gave us only implicit conditions and with so many general classes of integrable equations 
it is natural to try and work ones way in the other direction, i.e. reconstructing the integrable potentials. To 
demonstrate the idea let us first consider how to generalise the basic form of the equation in question. 

We want our main equation to assume some particular, known form with rational r(z), and notice immediately 
that energy enters in such a way that the infinity must be an irregular singular point. This could in principle be 
remedied by other parameters of y, but only for one fixed value of -B^, while the physical nature of the problem allows 
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any E to begin with (and usually and with a restricted, but still infinite spectrum). It is thus better to assume that 
infinity is irregular, and the special values of parameters for which it is not will up in the analysis anyway. 

At first one might ask whether there are rational potentials with the property that only the infinity is a singular 
point - like in the polynomial case. V{z) would then need at least one pole, not to be polynomial, and it could only 
be of first order and residue 1 if + V is to have no poles. The general potential would then be 



^ 1 

^ Z - Ck 



with a polynomial R{z), and r{z) reads 



k=l 



N , N 



r{z) = + R' + 2Ry —^ + 2y --E''. (28) 

^ Z-Cfc ^j-^ {x-Ck){x-ci) 

The only possibility of removing the poles is cancellation between the two sums, and expanding them at each Ck the 
condition is 



^ 1 

R{ck) + V 0. (29) 

^ Ck- ci 



1=1 



The poles are all different (otherwise the residue would not be 1), so that this system of N equations determines a 
polynomial of degree — 1 uniquely, but any other polynomial which vanishes at the poles can be added, so the most 
general form of R is 



R{z)=R^{z) + R-,{z)^{z-Ck). (30) 

k=\ 

where R\ satisfies (l29l) and i?2 is arbitrary. The system determining coefficients of R\ has no explicit general solution 
(inverse of the Vandermond matrix) for arbitrary iV, so let us again look at the potential with just one pole at zero 
so that 

n 

R{z) = zR^{z) ^ zY^Tkz'' , (31) 

fc=0 

and 

r{z) = 3i?2 + zi?2 + 2^^2 - (32) 

With no poles in r(z) and irregular infinity, only the first case of Kovacic algorithm can hold, so that the solution 
must be of the form (1111) . Substituting this into the differential equation and comparing with r(z) at infinity tells us 
that uj — zi?2 and also that the degree of P must be 1. Taking the polynomial to be monic P{z) = z +pO results in 
the final equation 

(£;2-2r2)p0 + i;2z = 0, (33) 

which can hold only when Ri is of degree 1 with po = ''o/''i and i?^ = 2ro (r\ = leads to = again). The 
integrable potential is then 

F(z) - i + z(ro + riz), (34) 
z 

and the solutions 

^\{z) = (riz + ro)exp (z^(3ro + 2riz)/6) , 

For higher N it is possible to proceed halfway in that the expansion of r{z) at infinity implies that lo — R{z), and 
the solution of the prescribed form P exp J uj leads to 

P" + 2RP' + (P? + R' - r)P = 0. (36) 

The highest order term vanishes when degP — N, so this is just a linear system for the coefficients of P{z), but to 
write it out explicitly one would have to know which amounts to solving (I29|) separately for each value of N. 
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V. 



Finally, let us look at a known type of equation for which nontrivial integrable cases exist. In contrast to the 
previous example, we expect r{z) to have a pole but want to restrict the degree at infinity (the polynomial part). As 
mentioned before, there is a special case of Heun equation which reduces to the confluent hypergeometric one and has 
a finite pole of order two in r{z), while only the energy accounts for infinity being irregular. Its general form is the 
Whittaker equation: 

and since its solvability conditions arc known, we can find out which potentials allow Liouvillian solutions. 
In order for the Dirac equation to assume such form, the potential has to satisfy 

V^^V..e'=(l-l,i!^. (38) 

which is a Riccati equation. Its rational solution can be obtained by the analysis of the poles of the right-hand side. 
The general solution to the above is 



V = m+ - 



d 

(39) 



The poles Cfc have to be included even though they do not correspond to the poles of V"^ + V because each terms 
of the sum is a solution to V"^ + V = 0. Since the pole at zero is essential in the equation, its residue is a 7^ 1. 
Alternatively, the above could be rewritten as 

V = m+- + ^, (40) 
z P(zj 

where P is a polynomial with d distinct roots Cfe, such that P(0) 7^ 0, and can be taken to be monic. 
Substituting this form of V to the coefficient r, one gets 

r = m^-E^ + ^(^ + g^+^(^^ + ")f + (41) 

2;^ z zP 

The restrictions on the singular points require that the last term is of the form a + b/z + c/z"^, so that we have an 
equation for P 

zP" + 2{a + mz)P' = {az + b+-)P. (42) 

z 

As the solution is supposed to be a polynomial, the coefficient a must be zero, or the terms of the highest order would 
not cancel. Similarly, the lowest order terms require that P{0)c = and since we excluded zero from the list of the 
roots c also must vanish. We now have 

zP" + 2{a + mz)P' = bP, (43) 
which is the associated Laguerre equation, with polynomial solutions only for 

b = 2md, dGN, m ^ 0, (44) 

so that 

P = Lf'-^{-2mz). (45) 
With such condition, the Dirac equation becomes 
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and assuming ^ E'^ (the special case will be dealt with below) a change of the dependent variable to 

u = 2\Jm^ -E'^z (47) 

gives 



, /I «(« - 1) TO a + , 
^ " = 7 + 2 + / o V', 



(48) 



so that the appropriate Whittaker's equation's coefficients are 

= -("- \) 

m (^^) 
K = =(a + a), 

and we simply have 

i'iz) = (2Vm2 ~ ^2^) . (50) 

Since the equation is of the second order, there will be also another, linearly independent solution given by 
W-n^^{—z). However, in application one will choose only one taking into account the physical context such as square 
integrability, and the domain of z. For example, in the z 3> region, one would not take the function W-k..^{—z) 
since it behaves as exp(z) asymptotically, while Wk,^{z) behaves as exp(— z). 

Here, however, we are interested in the general nature of the involved functions, rather than particular physical 
model, so we wish to impose other restrictions, namely the Liouvillian character of the solutions. For the Whittaker 
equation the conditions for that are well known Q and here they give 



which, upon substitution, leads to 



p /i + K — — G Z, or 
q:=/i-K-ieZ, 



1 = — (rv + d) + a £ Z, or 

Vm2 - E^ ^ ' 



(51) 



:(q! + d) - a € Z. 



(52) 



%/ m2 — E"^ 

In other words, when the potential is such that a — —d the solutions are always Liouvillian, and otherwise the 
parameters must satisfy 

lez, (53) 



Vm2 - E^ d + a 
which in turn means that the integrable part of the spectrum is bounded by 

£:2 , /d + a^^ 



0< — = 1- y— <1. (54) 

Note that when a is an integer itself, the case of / = ±a is excluded by to 7^ 0. 

For the particular case of to2 = £'2, equation is effectively the Bessel equation 

wV'H + wip'iw) + {w^ - (1 - 2af)ip{w) = 0, (55) 

when the variables are changed according to 

-0(z) = <p(w), w'^ := —8m{a + d)z, (56) 

so, assuming a ^ d, the solutions are given in term of the Bessel functions 



^{z) = v^Ji„2a(2v/-2TO(a + d)z), (57) 

and these are known to be Liouvillian only if their order, i.e. 1 — 2q; is half an odd integer. 
Finally, for a = —d, the equation is solvable by elementary functions of 

i^(z) = ciz" + C2zi-". (58) 
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VI. 



To put the results together, we have shown that the one-dimensional Dirac equation with a rational potential will 
lead to integrability conditions which in general will be given by polynomials defined recursively. Also, that rational 
potentials can be determined if some natural conditions are required of the resulting equation itself. The direct and 
inverse results are as follows. 

The simplest form ofV = a/z + m + Xz gives rise to the biconfluent Heun equation (if none of the parameters is 
zero). Its solvability is determined by the polynomial condition p^i{a, m) = 0, obtained by one of the relations (|16|) 
(together with the respective conditions on E, X and d given by ([H]) ) - if any of the four can be solved it gives a 
Liouvillian solution. However, we are only interested in real values of the parameters and the three theorems [1] [2] and 
[3] characterize (to some extent) the position and asymptotics of the real roots. 

The special cases of the above potential, for vanishing parameters are 

1. a = 0, the Dirac oscillator solvable by Hermitc functions @. 

2. TO = 0, the Laguerre equation, solvable with Ln ^^^{\/—Xz), with 4An = —E^. Liouvillian when n or n + a— 1/2 
is an integer. 

3. A = 0, the Whittaker equation of section V, summarized below. 

If we restrict the problem somewhat by requiring that r(z) remains a polynomial, although V{z) might still be 
rational, we obtain a general condition that the potential is of the form 

N , N 



Viz)^y^^— + Ri{z) + R2{z)Y]{z~Ck), (59) 

with i?i of degree — 1 such that 



Ch 

k=l " k=l 



^i(cfc) + V =0, (60) 

l^k 

and i?2 and arbitrary polynomial. This is the most general rational potential that leads to a polynomial coefficient 
r{z). Due to Ri being given only implicitly (although uniquely) only the one-pole case in analysed in section IV. The 
potential turns out to be 

V{z) = ^ + roz + riz'^, (61) 

and the integrable values of the energy are E^ = 2ro- 

Finally, choosing a predetermined r{z) which is simpler than the Heun equation, but still has infinity as an irregular 
singular point and one other pole, leads to the Whittaker (confluent hypergeometric) equation, and the potential to 
be 

a P'iz) 
z F{z) 

where P is expressed by the associated Laguerre polynomial of degree d as 

P(z) = Ll''-^{~-2mz), (63) 

and m is necessarily non-zero. The solutions are then given by the Whittaker functions Wk.^^{z), and are additionally 
Liouvillian (can be expressed in closed form) when 

1. a = —d regardless of energy, or 

2. 2a e i + Z for E"^ = m'^, 

3. for particular levels satisfying -^=^=5= (a + d) ± a £ Z. 



12 



a, E 


F, <= r 




^2 ^ ^2 


a G C 




W-«,„(-u) [L] 


\/iJi-2c«(w), yiYi_2ti(iy) 


2a - i G Z 




[L] 


V^-^l-Zalw), yiyi-2a(w) [L] 


a = —d 


[L] 







Table I: The list of bases of solutions depending on parameters. W denotes the Whittaker function, J and Y the Bessel 
functions. The auxiliary variables are: w = 2-^— 2m(a + d)z, u = 2\Jm?- — E^z. [L] signifies a Liouvillian pair of solutions. 



The third possibility gives countably infinitely many energetic levels which are confined to _B G [— ™, ni\. All the cases 
are put together in table ID 

It should be clear from the general expression of the potential (|62| that it describes a very wide class, with an 
arbitrary number of singularities, thus modelling drastically different physical situations. The mathematical analysis 
given here gives the full information about the integrability and formal solutions, but it is beyond its scope to select 
the appropriate eigenstates and respective energies corresponding to a particular quantum system. However, just as is 
the case with the Dirac oscillator, the integrability condition can be thought of as some quantization rule. As shown 
in dual rules will hold for the vector coupling upon exchanging the roles of E and m. Either way, once can speak 
of some form of quasi-integrability here, since Liouvillian solutions depends on values of (at least) two independent 
physical parameters, so even before taking the physical restrictions into account, we are left with at most countably 
many explicit solutions. 

Concerning the property of shape invariance, the recovered potential contains the known Coulomb case (albeit the 
spatial variable is then radial) listed in Q. This is the only case when the property holds, so that the whole class is 
in fact new, and the method described can be used to uncover more solvable cases by loosening some of the initial 
requirements. Future generalisations might include target equations with more singular (regular or not) points, or 
simply integrability conditions for given, more complex potentials (trigonometric, meromorphic etc.). 

Of course the Liouvillian solutions are not the only notion of integrability, and most linear equations can be solved 
by means of series or numerically, but the concept of an exact solution remains as the starting point in any such 
approach. Any further study of the properties of the system in question, especially qualitative, it greatly facilitated 
by working on explicit formulae, which this work hopefully provides. 



Appendix A 

When the degree d = 2n is even, p_i will be an odd polynomial, so m = will make it trivially zero. Also, writing 
the relation for even I and taking the limit we get (in the first subcase) 

lil + 2a-l)pi^~2{d~l + 2)pi^2, (Al) 

so that only even indexed polynomials are generated from pd — I, leading to po but not p-i- It is thus necessary to 
consider polynomials r := pi/m, which changes the relation to 

l{l + 2a- l)n + 2{l ~ l)pz_i = -2(d - I + 2)r,_2, (A2) 

Since m does not enter the relation the limit of to changes only the initial polynomial r^-i = —d. Note that r; 
is defined only for odd / and is the same as above, with its index even. The energy forces d to be at least 2, and 
then r_i = 2a/3, so there is one simple root (albeit excluded). 

The d > 4 case can be dealt with inductively. First the explicit formula for the even indexed polynomials is 

-2fc(« + fc-l/2) 

k=l+l/2 

So that pd-2 = ~d{2a + d— l)/4 has only one root, and as / decreases there appear new simple roots at half (negative) 
integers, given by 

ak = ^-k, fc = ^,...,l + ^, (A4) 

so that pi all share the first zero, and each has {d—l)/2 zeroes altogether. The sign of the polynomial does not change 
past its final zero, and can be determined from (IA3|) 

sgn(p,(a,,M + l/2)) = (-l)('^-')/2, (A5) 
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where the subscript in ai^j^i indicates that the maximal root depends on /, and the reason for taking the value 
specifically at ai^M + 1/2 will become clear in a moment. 
The initial inductive step will be ^ = d — 1 with the relation 

(d - l)(d + 2a - 2)(-d) + 2(d - 2)pd-2 - -6rrf_3, (A6) 

where r^^i = —d, and the other polynomials are linear in a but it is more instructive not to substitute them explicitly. 
Now, at the first (and only) root of Pd-2, Oio^ the above simply reduces to 

1) = -6rd-3(ao) > 0, (A7) 

whereas at ao + ^ we have {d + 2ao — 2) = so 

2{d - 2)pd-2{ao + 1/2) = -6rd_2(ao + 1/2). (A8) 

By (jA5|) . the above translates to 

(_^)(d-d+2)/2 ^ _^ ^ -sgn(rd_2(ao + 1/2)) = -sgn(rrf_2(ao + 1)) < 0, (A9) 

where the last equality follows from the fact, that rd-2 is linear (by construction) and already changes sign once in 
the interval (ao,ao ~^ \)- Consequently there is a root, and since ao + 1 is the next zero in the sequence (jA4p . we 
know the signs of rd-'i at the first two ak- 

The next step is to show that all following n have the same signs at the (d — I + l)/2 roots of in alternating 
sequence, and that there is an additional change of sign at the (d— ^)/2-th root of p;_3. This is to extend the previous 
paragraph, where the signs agreed on the first root, and there was a change at the next one. For / < d — 1 the relation 
at any zero Uk of p;-i is 

l{l - 2k)ri(ak) = -2(d - I + 2)r,_2(afc), (AlO) 

so that the signs agree because 2 + I < 2k < d. Thus as the signs alternate for r/ they also alternate for r/_2 
indicating (d — ^ + l)/2 simple zeroes. Note that at the first zero ao, the signs of all r/ are thus —1 and accordingly 
sgn(n_2(a;_i_M)) = (— l)(''~'+i)/2. As before we look what happens past the maximal zero at a/_i_M + l/2 — {l — l)/2 
so that the first term of the relation is zero and 

2{l - l)p,_i(az_i,M) = -2(d - I + 2)ri_2{ai-i,M), (All) 

or, in term of signs, 

(_l)(d-i+i)/2 ^ _.sgn(r,_2(a/-i,M + 1/2)). (A12) 

n_2 is at most of degree {d — I + l)/2 as a sum of polynomials of such degree, and the last change of sign means this 
is strictly the degree with all roots real and simple. With no more sign changes it also follows that r/„2 («;-!, m + 1) = 
(_l)(rf-'+i)/2 and since that is the last root of p;_3 the last induction step is completed. 

The last polynomial, r_i, must thus have d/2 real simple zeroes, although the last one is a = 0, for, as mentioned 
before, p_i is proportional to api. 

The third subcase is very similar to the first one, with positive A this time, and the even pi characterized by 

^ fc(2a-2fc-l) 

P'= 11 d-2fc + 2 ' 
fe=i+;/2 

ofe = (d+l)/2,...,(/ + 3)/2, (A14) 
sgnipiiai,m - 1/2)) = (-l)('^-')/2, (A15) 

where ai^m is the smallest root and the greatest root is shared by all pi. For odd I the relation reads 

l{l + l- 2a)ri - 2lpi^i = 2{l - d - 2)n_2, (A16) 

and exactly as above the signs of ri and r;_2 agree on zeroes of pi-i. The smallest root of is 1 + //2 so the 
relation at ai-i^m — 1/2 gives 

Ipi^iiai^um - 1/2) ^{d-l + 2)n_2(a(-i,m - 1/2). (A17) 
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We can then proceed inductively, as before, by noting that r^-i = d so that all ri are positive at the maximal root, 
and they alternate signs at a^, with the additional change of sign resulting from the above equation. The only trivial 
case to discard is d = for which r_i — 1, while for d > 2 it follows that there exist d/2 real simple roots. 

Subcase 2 differs in two ways. Firstly, the sign of A is arbitrary, but this only means that for A < the signs 
of subsequent r; will be opposite at zeroes of pi instead of the same. That will not change the fact that the signs 
alternate and account for all but one zero. Here is where the second difference appears, for we have 

fc(2a + 2fc-l) , 

k=l+l/2 

afc = (l-d)/2,...,-(/ + l)/2, (A19) 

and as the relation is now 

l{l + 2a- l)ri ~ 2{l + 2a - l)p;_i = 2A(; -d- 2)r,_2, (A20) 
it means that at ai-i^M + 1/2 we get 

= 2A(/-d-2)n_2((l-0/2)- (A21) 

So that the last zero of ri is known explicitly and it ensures the desired changes of signs at roots of the appropriate 
pi. Since r^-i = 2a + d and pd-2 = —d{2a + d — l)/2 have different roots (and hence the sign changes accordingly), 
induction on I yields the existence oi l + d/2 roots of r_i with the explicit root being a = 0. Note that this also holds 
for d = where r_i = 2a. 

Finally, in subcase 4, there is a complication due to r^-i = —2a + d + 1, pd-2 = d{—2a + d + l)/4, so that all 
polynomials share a root. As mentioned before, it corresponds to the degenerate E — case. This is the greatest of 
the roots of pi which are 

k(-2a + 2k + l) , , . 

n ^^(Z32fcT2r'"^-(' + 3)/2,...,(d+l)/2. (A22) 

Let us consider A > for clarity (the negative case is completely analogous). Note that pi will be positive for a less 
than the smallest zero because the leading coefficient will be (— a)^''"'-'/^ and since the relation is 

+ 2a)ri + 2{l - 2a)p;_i ^ 2{d - I + 2)r;_2, (A23) 

the same will be true for r; and its smallest zero. We will see that apart from a = (c? + l)/2 all the zeroes of r;_2 lie 
between the zeroes of pi^'i- Direct computation shows that the smaller zero of r^-a is ?,d{d — l)/(6d — 4) which lies 
between {d— l)/2 and (d+ l)/2. In particular r^-a is positive at {d— l)/2 — ad-i.m- Assume that this is true for r; 
and the relation for zeroes of pi-i other than {d + l)/2 gives 

l{l + l~2ak)ri{ak)^2{d~l + 2)ri^2{ak), (A24) 
with the coefficient of n being negative due to (|A22[) . At ai^i ^ — 1/2 = (Z + l)/2 the relation reads 

- 2p,_i(ai„i^,„ - 1/2) = 2{d-l + 2)n_2(a;-i,m - 1/2) < 0, (A25) 

but at a;_3_m — 1/2 we get 

/n(a,_i,™ - 1) = 2(d - / + 2)rKa;-i,™ - 1) > 0, (A26) 
because r; is positive to the left of its smallest zero, wliicli is grcciter tlis-n (y.i—\_^rn- 

In other words, ri-2 is positive 

at a;_3^m, and negative at a;_i^m and then changes the sign between the subsequent zeros of pi-i because n does 
so. At the last but one zero (of pi-i) it thus has the sign of (— l)('^~'~i)/2 (number of zeroes of pi-i except the last), 
and close to the left of the last zero pi-i and ri both have the sign (— l)(''^'^i)/2 because their leading terms are 
(_Q,)('i-'+i)/2_ However, they enter the relation (jA23p with negative coefficients (/ + 1 — 2q; < when a > d/2) which 
means that n_2 changes sign additionally before the last root of {d+ l)/2. Last thing to check is the d = case, for 
which r_i = 1 — 2a. This finishes the proof that r_i has 1 + d/2 real simple roots. 
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